Diffusion of small particles in a solid polymeric medium. 
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We analyze diffusion of small particles in a solid polymeric medium taking into account a short 
range particle-polymer interaction. The system is modeled by a particle diffusion on a ternary lattice 
where the sites occupied by polymer segments are blocked, the ones forming the hull of the chains 
correspond to the places at which the interaction takes place, and the rest are voids, in which the 
diffusion is free. In the absence of interaction the diffusion coefficient shows only a weak dependence 
on the polymer chain length and its behavior strongly resembles usual site percolation. In presence 
of interactions the diffusion coefficient (and especially its temperature dependence) shows a non- 
trivial behavior depending on the sign of interaction and on whether the voids and the hulls of 
the chains percolate or not. The temperature dependence may be Arrhenius-like or strongly non- 
Arrhenius, depending on parameters. The analytical results obtained within the effective medium 
approximation are in qualitative agreement with those of Monte Carlo simulations. 

PACS numbers: 05.60. Cd; 82.35.Lr; 05.45.-a 



INTRODUCTION 

The literature treating the problem of diffusion of small 
molecules in solid polymeric media is surprisingly lim- 
ited when compared with the huge amount of results ob- 
tained for diffusion in solids in general. This is even more 
surprising if one takes into account the enormous role 
polymeric materials play as encapsulants and isolating 
materials in technical devices. Furthermore, the ther- 
modynamics of polymeric solutions and the dynamics of 
polymers in solutions are well understood since the sem- 
inal works by Flory [3 and Huggins [4 . The situation 
usually considered is the one in which polymer molecules 
constitute the solute of the solution. Only few times the 
roles have been inverted and polymers have been taken as 
solvent molecules forming a matrix in which small solute 
particles are let to diffuse. Early works done in this di- 
rection [1, 2 analyze the concentration dependence of the 
small solute diffusivity through experimental adsorption 
and desorption curves. In particular, in [2j, Fujita con- 
cludes with the necessity of theoretical and experimental 
investigations of the "characteristic differences" between 
the cases of a good or bad solubility. Up to the authors' 
knowledge, this request has not been satisfied yet. Aim 
of the present work is to give a partial answer by provid- 
ing a qualitative analysis of the way the particle-phobic 
or particle-philic nature of the polymer chains affects the 
diffusion process. 

More specifically, we consider a set of particles dif- 
fusing in an amorphous solid polymeric medium, in a 
model being a close relative of a classical Flory-Huggins 
model of polymer solutions. In the present work we 
adopt the ternary lattice representation corresponding to 
a polymer-solvent-void system close to the one proposed 
in Ref.[6 : In the two variants of the model considered 
we take a site of a lattice to represent a polymer seg- 
ment, an interaction site in the vicinity of a segment, or 



to be empty. The concentration of solute molecules is 
considered low, and their interaction with each other is 
neglected. 

In the first variant of the model polymers are rep- 
resented by chains of occupied sites and their nearest 
neighbors are considered as interaction sites. Sites not 
belonging to either of these two categories are considered 
as voids. This lattice model is exactly the one we use 
in simulations. Analytical calculations refer to a sim- 
pler mean field Flory-Huggins-like model, built by dis- 
assembling the chains and letting polymer segments, in- 
teraction and empty sites fill the space in a completely 
random fashion at given concentrations. The situations 
are discussed in depth in Section 2. Details of analyti- 
cal calculations are given in Section 3 with a particular 
attention to the variations to be made with respect to 
the conventional effective medium technique. In Section 
4 the interaction between the polymers and the small so- 
lute molecules is temporarily switched off and the model 
is reduced to a pure percolation problem in the presence 
of polymer chains. This is done in order to estimate the 
error introduced by the mean field approach and the de- 
pendence of the diffusion coefficient on the chain length. 
Section 5 is devoted to the role of interaction sites, and 
Sections 6 contains our conclusions. 

THE MODEL 

We model our solid polymeric matrix by a three- 
dimensional cubic lattice on which the chains are mod- 
eled as phantom random walk chains of length This 
chain conformation corresponds to the Gaussian nature 
of chains in melts from which our solid matrix is obtained 
by quenching. The whole matrix is considered as static: 
no chain motion is taken into account. After the system 
is created, the sites of the lattice occupied by chains are 
considered impenetrable for small solute molecules. The 
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whole system is then modeled by a ternary random po- 
tential landscape. The corresponding lattice is outlined 
in Figure [l] 

The sites occupied by polymer segments are impene- 
trable for solute molecules (hard core interaction, inter- 
action energy U = 00) and represented as black sites 
in Fig. [1] The number concentration of these sites is 
^3 = M3/M where M3 is their total number and M is 
the volume (total number of sites) of the lattice. 

The particle-polymer interaction is considered to take 
place only if the molecule occupies a site which is a 
nearest neighbor of the one occupied by a polymer seg- 
ment. The particle-polymer interaction at these sites cor- 
responds to the interaction energy U = e whose sign fixes 
the nature of the force experienced by the particles: if e is 
negative, this interaction is attractive; if e is positive, the 
interaction is repulsive. These interaction sites are repre- 
sented in red in Figure [l] and their number concentration 
is ^2- 

Remaining sites are considered as simple voids with en- 
ergy U = and where particles perform a free motion not 
being subjected to any force. The number concentration 
of these sites is ^1 = 1 — ^2 — ^3 and they are represented 
in white. Our system is thus represented by a random 
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thermodynamical properties of polymeric solutions, in 
which the number concentrations of the sites occupied 
by polymer segments is kept, but the correlations be- 
tween their positions (necessarily introduced by the ex- 
istence of chains) are fully neglected. This model corre- 
sponds to filling the lattice at random with black, red and 
white sites at given number concentrations. In this way 
each lattice site is assigned an energy value Ui which can 
take one of the three values 0, £, or 00 at random, with 
probabilities The existence of an infinite cluster of 
black sites, which we need to preserve the solidness of 
the system, is guaranteed by taking the concentration 
(j)s above the percolation threshold which is known to be 
approximately 0.32 for the three-dimensional simple cu- 
bic lattice. We denote this construction as mean field 
lattice and represent it in Figure [2] The diffusion on 
this mean field lattice is then treated using the effective 
medium approximation for a diffusion in a random po- 
tential landscape, as discussed in Section 3. The mean 
field / effective medium results are compared with the 
results of direct numerical simulations discussed above, 
and show qualitatively similar behavior. 
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FIG. 1: (color online) The ternary lattice of the initial model 
and different kinds of cells with corresponding concentrations 
and energy values. 



FIG. 2: (color online) Mean field ternary lattice and different 
kinds of cells with corresponding concentrations and energy 
values. 



(but correlated) ternary lattice with the sites assigned en- 
ergies Ui which take the values 0, e, or 00 for the white, 
red and black sites respectively. In this medium the small 
molecule diffusion is numerically simulated as a nearest- 
neighbor random walk with transition rates between the 
sites given by the corresponding energy differences: 



The constant rate wq defining the time unit of the process 
is set to unity in all simulations, P is the usual 1/KbT 
term and Kb is the Boltzmann constant. 

The analytical calculations are performed within a 
simplified model which strongly resembles the classi- 
cal Flory-Huggins model (^16^) used for description of 



EFFECTIVE MEDIUM APPROXIMATION FOR 
DIFFUSIVITY 

The particles' motion in a random potential landscape 
is described via the usual master equation 

Qi = ^{^ij(lj -WjiQi), (2) 
j 

where Qi is the probability for a particle to be at a site 
i at time t and Wij is the transition rate from site j to 
site i given by equation ([T]) for j and i nearest neighbors 
and equal to zero otherwise. For the sake of generality 
calculations will be referring to the (i- dimensional case. 

We multiply both sides of equation ([2| by the number 
of particles N and obtain the master equation for the site 
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mean number or "concentration" function rii = Nqi. 

ni = ^ {wijUj - WjiUi) . (3) 

3 

Assuming the existence of an equilibrium state, the tran- 
sition rates are naturahy hnked through the detailed 
balance condition at equilibrium Wijn^ = wjiU^^ where 
= Nq^ and ex exp {—/3Ui) is the equilibrium proba- 
bility to find the particle at site i. Thus one can introduce 
the symmetrized rates gij being the properties of a bond 
of a lattice, 



with 
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(5) 



where ^((/), e) is the normalization factor of q^ (the parti- 
tion function for the small particles equilibrium distribu- 
tion) and (j) is the triplet (^1,^25^3)- Then the analogy 
between the diffusion and the electric conduction in a 
random medium can be used [71 [3 [M] : the correspond- 
ing diffusion coefficient is connected with the macroscopic 
conductivity {g)em of a disordered lattice with bond con- 
ductivities gij via ^14 | il7 j 
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with a the lattice spacing. Our system exhibits four dif- 
ferent bond conductivity values depending on the color 
of the sites involved. These are 

gi = go; g2 = goe~^'; gs = goe~^'^'^; g^ = 0. (7) 

Figure [3] gives an overall view of this situation. 
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FIG. 3: (color online) Bond conductivities for the correspond- 
ing site couples. 



The effective conductivity {g)em can then be calcu- 
lated within the effective medium approximation (EMA). 
There is however a subtlety in application of the effective 



medium approximation to site models like ours. The gen- 
uine continuous EMA of Ref. [11 and its lattice variant 
Ref. [10 describe well the behavior of bond percolation 
model but fail to reproduce the behavior for site perco- 
lation to which our ternary model reduces when e = 0. 
For the site model EMA procedures to obtain gem were 
proposed by Bernasconi and Wiesman in ^ and Yuge 
in [16 : in both the usual effective medium procedure is 
varied in order to take into account the correlations be- 
tween consecutive bonds naturally arising in this kind of 
systems. The fact that such correlations arise is easily 
understandable when considering a simple example: let 
us take three neighboring sites j and k where i and k 
are two different nearest neighbors of j, and consider the 
bonds ij and jk. If j is, say, white none of these bonds 
can have conductivity ^2 because both of them involve 
the white site j, and their conductivities are not inde- 
pendent as assumed in the bond-based approach, and the 
usual EMA technique has to be appropriately changed. 
Our approach here follows the lines of Ref. [16 . 

Calculations start by considering for any site of the 
lattice possessing a color index a = 1,2,3 (corresponding 
to white, red and black, respectively) the mathematical 
expectation of the conductivity g^ of a bond starting from 
it: 



gi = (Pigi + 02^3 

g2 = hgs + ^2^2 
^3 = 



(8) 



These values appear in the system according to the the 
probabilities of their respective sites 



(9) 



i=l 



The effective conductivity is then obtained through the 
usual self-consistency condition [lOj 



{d - l)gem + g 



0. 



(10) 



where d is the dimension and {■) p is the average with 
respect to the distribution P above. If we now define a 
rescaled effective conductivity 
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and introduce the arithmetic mean and the 0/e-weighted 
average of the quantity Ei = e~^^^/^ 



E= -(1 
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) and (E) = ^1 + (j)2e 
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equation (10) reduces to a quadratic equation for /em, 
/L + /em + c(<?,e) = (13) 
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with 



b{$,s) = {E)(dE-d{E)) 



c(0,e 



The value of Dem follows from the solution of this equa- 
tion via 
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where I^o = cl'^wq is the diffusivity of a lattice where 
all sites are white. The critical threshold at which 
vanishes can be obtained by setting c(^, e) = 



1 - 1/d. 



(15) 



At variance with a classical binary Flory-Huggins sit- 
uation, even in a random model with given ^3 it is hard 
to get an analytical estimate for ^1 and ^2 as functions 
of the known ^3 concentration. In a random case 02 
is the total perimeter density of black clusters in a site 
percolation (whose behavior is in principle known, but 
whose values have to be estimated numerically), and (j)i 
corresponds to the rest of the sites. This task gets even 
harder if the chains are present, and the numerical sim- 
ulations show that the existence of the chains does mat- 
ter. To overcome the problem, we simulate our polymer 
model first and extract the numerical values of 0i and 
02 from these simulations. These numerical values are 
then used in the corresponding EMA calculations, whose 
predictions, in their turn, are compared with the results 
of simulations of diffusion. 



PURE PERCOLATION (BINARY) MODEL 

It would be nice to know, how large is the typical er- 
ror arising from disregarding the chain structure of black 
sites, and what is the role the chain length plays in the 
simplest case, namely in a percolation model with corre- 
lated black sites given by the chains. In this model the 
red and the white sites are indistinguishable, they have 
the total number concentration 01 = 1 — 03, and the 
result of our previous consideration reduces to the origi- 
nal Yuge's result for site percolation. This is exactly the 
situation discussed in the present section. 

Thus we consider a pure percolation situation in which 
the only interactions are the excluded volume ones and 
our lattice consists of only black and white sites, red ones 
are absent. The results of simulations for the systems of 
chains of different lengths are shown in Fig. |4ja). The 
figure representing the dependence of the diffusion coeffi- 
cient on the concentration of sites occupied by segments 
of the chain shows this for the chain lengths from / = 1 



(usual Bernoulli site percolation problem) to / = 10. The 
simulations were performed also for longer chains, but 
for / larger than 10 the corresponding graphs are indis- 
tinguishable from that for / = 10 within the statistical 
accuracy. Thus, a result for / = 100, (not shown) is in- 
distinguishable from the one for / = 10 on the scales of 
Figga). 

Details about simulations are readily given: Simple 
random walks of / steps are let run independently in a 
lattice of 400^ sites with periodic boundary conditions. 
This operation is stopped when the total segment concen- 
tration of segments (sites visited at least once) is within 
0.01 from the desired value of 03. Once the environment 
is created, 10^ random walks of 10^ to 10^ steps, depend- 
ing on the speed of homogenization of the system, are 
launched from a free site chosen at random in a cube of 
50^ sites placed in the center of the medium. With this 
choice, the probability for a diffusing particle to reach 
the borders of the lattice is extremely low and doesn't 
spoil the statistics. The algorithm used is the Monte 
Carlo Blind Ant one. The whole procedure is then re- 
peated for 10 different lattice realizations and averages 
are taken. We have observed a normal diffusion process 
(r^(t)) (X t from which the proportionality constant Ds 
has been extracted and reported in Fig. [4]^a). The ho- 
mogenization of {r'^it)) slows down in the proximity of 
the critical point. For this reason 10^ time steps become 
insufficient and the diffusivity is systematically overesti- 
mated. Our attention however is focused on a range of 
values of 03 which are above the percolation threshold. 

The curves do not differ drastically, but definitely show 
different percolation thresholds 0i(O depending on I. For 
the Bernoulli case the total behavior of diffusivity is re- 
produced sufficiently well by EMA for 03 close to unity 
but departures from the EMA line for concentrations 
close to a critical one. For longer chains the critical 
concentration gets lower, and the diffusion coefficient at 
given 01 gets larger than for the Bernoulli case. Although 
different, the curves however show a large amount of uni- 
versality which is unveiled when rescaling the concentra- 
tion and diffusivity according to 

0; = |-1 and 4 = D3^i-M, (16) 

SO that the critical concentration is mapped onto the 
point (j)^ = 0, seeji^b). In this case all the curves fall 
onto the same master curve, and the mean- field result, 
rescaled accordingly, gives a straight line (of slope 1) 
which reproduces the results of simulations astonishingly 
well up to the critical domain. This high degree of univer- 
sality shows that the correlations introduced by the ex- 
istence of the chain are not of high importance and can 
be fully accounted for by rescaling the results of EMA 
according to the equations above. The corresponding 
critical concentration has however to be obtained numeri- 
cally. Alternatively, it can be extrapolated from the slope 



5 



CO 

*Q 0-4 



- ; = 1 

■1 = 2 
■1 = 6 
■1 = 10 

- EM A 




1.0 -0.5 0.0 0.5 



1.5 2.0 2.5 



(a) 



(b) 



FIG. 4: (color online) (a) Normalized diffusivities in the 
purely percolation case vs. the concentration 01 of white sites; 
(b) Rescaled normalized diffusivities vs. rescaled white con- 
centration 01. The dotted line represents in both figures the 
effective medium approximation. 



of diffusion coefficient for concentrations close to unity. 



RESULTS FOR TERNARY MODEL 

In this section we discuss results for the normalized 
effective diffusivity Ds{(j),e) and concentrate on the role 
of interaction energy £ between the diffusing particles 
and the polymer matrix. All the figures refer to the 
three-dimensional case. The reduced interaction energy 
/3e = e is chosen to span in the interval [—5, 5] accord- 
ing to the following reasoning: typical absolute values of 
^xx/Kb^ the coupling strength of a Lennard- Jones po- 
tential describing the interaction between two atoms of 
the same kind X, can be roughly enclosed in the inter- 
val corresponding to temperatures [0, 500i^]. In order to 
consider the interaction between two different atoms X 
and Y the Lorentz-Berthelot mixing rule is used to ob- 
tain sxY = ^/sxx^YY which, being an average, belongs 
to the same interval. Using e in place of £xYi considering 
both positive and negative values and taking the temper- 
ature not too far from the ambient one, it is straightfor- 
ward to see that the choice £ G [—5, 5] is a reasonable one. 
For the discussion of the Arrhenius-like or non-Arrhenius 
temperature dependencies in Sec. [] broader bounds are 
used, ee [-10,10] 



Effective diffusivity vs interaction energy 

Let us first discuss general features of the dependence 
of the diffusion coefficient on number concentrations and 
on interaction energy e. The EMA results for Ds{(t)^ e) for 
the three different cases corresponding to different rela- 
tions between and (j)2 for (ps fixed are shown in Fig. |5j 
These plots show the behavior for the attractive and re- 
pulsive interaction and the way the diffusivity approaches 



zero when the black sites concentration approaches its 
critical value 4>3 = 2/3 (see Eq.(15)). At this value in 
fact, particles remain confined in finite subregions of the 
system, due to the overwhelming predominance of poly- 
mer segments. 

Plots are given for three different sets of the (j)k values 
in order to consider symmetrically the situations in which 
red sites are in minority, equally probable or predominant 
with respect to the white ones, at given ^3. For this 
purpose we introduce a real parameter 7 G [0, 1] defining 
the number concentrations of white and red sites as 



h = 7(1 - ^3) 

h = (i-7)(i-^3). 



(17) 
(18) 



Graphs are then taken for three different values of 7 
(color online): 7 = 3/4 (blue dotted lines, (j)2 < (/>i), 
7 = 1/2 (black dashed lines, (j)2 = (/>i) and 7 = 1/4 
(red solid lines, ^2 > ^i)- This imbalance will deeply 
influence the behavior of the effective diffusivity when s 
crosses the zero value. 

In the symmetric case (j)2 = ^3 is invariant un- 
der the change of the sign of interaction energy e ^ —£. 
On the contrary, when the white-red balance is broken, 
the effective diffusivity decreases or increases depending 
on the sign of the energy parameter and on the value of 
7. Let us consider the situation in which (j)2 > 4>i (e.g. 
7 = 1/4, Fig. jsjc)) and restrict our attention on the 
attractive e < region; with this choice, we increase the 
number of the red-red ^2 bonds (showing larger conduc- 
tivity) with respect to the number of the white-white qq 
ones which have the lowest conductivity. This results in 
a global increasing of the effective diffusion constant. If 
we now invert the sign of £, i.e. consider the repulsive 
interaction, the ^2 bonds will still be the most numer- 
ous, but now have lowest conductivity value, decreasing 
in this way the whole diffusivity of the system. The op- 
posite happens if we consider (j)2 < 4>i] the corresponding 
graph in Figjsjb) is a mirror image of the one in Fig. [sjc). 
The comparison between the mean field calculations and 
the Monte Carlo simulations performed in the original 
ternary lattice corresponds to the chain length / = 100 
is quite satisfactory (figjsjd)). Once the desired poly- 
mer concentration (j)^ 0.4 was reached and the poly- 
mer matrix is set up, the energy value £ is assigned to 
all the nearest neighbors of the segments and their con- 
centration (j)2 is measured. All results are averaged over 
10 realizations of the polymeric matrix. In each of them 
random walks of 10^ steps were performed as described 
above. The number of realizations per energy landscape 
is 10^ times. The numerical values of number concentra- 
tions are ^3 = 0.39 and ^2 = 0.43. The numerical result 
is then plotted together with the mean field calculations 
in which the same values are used (figure [sjd)). We note 
that the value of the polymer concentration is close to 
the critical domain in Fig. |4ja) corresponding to e = 0, 
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FIG. 5: (color online) Mean field normalized effective difFu- 
sivity Ds vs e in the case 03 = 0; 0.2; 0.4 and: (a) 7 = 1/2; (b) 
7 = 3/4; (c) 7 = 1/4. (d) Comparison between theory and 
simulations (black dash-dotted line) in the case 03 = 0.39, 
02 = 0.43 (7 = 0.295). 



behavior can be explained as follows. Let us focus our 
attention again on the red (solid) line in the negative e 
half-plane. Under segment concentration 03 = 0.4 the 
black infinite cluster exists but is not dense enough to 
prevent the existence of infinite white or red ones. The 
concentration of red sites is 02 = 0.45 (7 = 1/4) and 
thus lays above the percolation threshold for a cubic lat- 
tice. This means that red sites form an infinite cluster 
crossing the whole system, and once a particle finds it, 
it can travel on it through the whole system rather than 
escape from it by activation. As a consequence diffusivity 
saturates and the system never freezes. In the repulsive 
region, the same behavior is shown by the blue (dotted) 
line, indicating the existence of a white infinite cluster. 
The black (dashed) line, the one for symmetric situation 
03 = 0.4, 02 = (/>i = 0.3, doesn't show any saturation. 
This suggests that in such a case white and red con- 
centrations are below the percolation threshold, and the 
activation processes are necessary to traverse the system. 

Figure [6]^b) shows the comparison between theory and 
simulation Arrhenius plots in the original interval e G 
[-5,5]. 



SO that the total accuracy of EMA is not too high in this 
domain. However, the EMA-results reproduce the de- 
pendence qualitatively well, and, moreover, the accuracy 
of EMA improves for higher interaction strengths. 



Arrhenius vs. non-Arrhenius behavior 

A non-trivial aspect of the dependence of diffusivity 
on the interaction strength is revealed by the Arrhenius 
plots shown in figure [G^a) where the logarithm of is 
plotted as a function of e = e/KbT in the wider inter- 
val [—10, 10], to investigate the role played by activation 
in the diffusion process; the segment concentration is set 
here to 03 = 0.4. The three curves in Fig. |6|a) corre- 
spond to the values of 7 = 1/4, 7 = 1/2 and 7 = 3/4. 
As in the previous figures, the curve for 7 = 1/2 repre- 
sents an even function of £, and the curves for 7 = 1/4 
and 7 = 3/4 are mirror images of each other. For e close 
to zero, the activation process is not relevant, the curves 
fall together and reproduce the diffusion constant in the 
black-and-white lattice of section 4. When moving away 
from the £ = value, the activation acquires importance. 
For 7 = 1/2 this behavior becomes Arrhenius-like and the 
curve shows a linear decay for both signs of e provided 
the interaction is strong enough. For asymmetric cases 
7^1/2 the Arrhenius behavior is seen only for interac- 
tion energy of the corresponding sign (attractive interac- 
tion for 7 < 1/2 and repulsive interaction for 7 > 1/2). 
For the opposite sign of interaction, at low temperatures, 
or high absolute values of £, the lines become horizon- 
tal, quitting the Arrhenius regime. This non-Arrhenius 




(b) 

FIG. 6: (color online) (a): Arrhenius plots of the different 
mean field normalized effective diffusivities vs. £ at 03 = 0.4. 
(b): Comparison between theory and simulation (black dash- 
dotted line) Arrhenius plots at 03 = 0.39, 02 = 0.43. 

On the total the following regimes of behavior can be 
qualitatively distinguished: 

(1) If the concentration of black sites is so high that 
percolation on red and white sites is not possible, the 
diffusion coefficient vanishes. 
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In the case when percolation over the red- and- white 
domains is possible, the diffusion coefficient is nonzero, 
and its behavior as a function of temperature depends on 
the percolation properties of red and white clusters, and 
on the sign of interaction energy. 

If the interaction is repulsive, two regimes appear: 

(2) If white sites percolate, the diffusion over the white 
cluster is always possible and does not need activation. 
The temperature dependence saturates. 

(3) If white clusters do not percolate, the diffusion is 
only possible over red sites, and involves an activation 
process; its temperature dependence shows the Arrhenius 
behavior. 

In the case of attractive interaction the roles of white 
and red sites interchange, and percolation over red sites 
is what determines the temperature dependence of the 
diffusion coefficient: 

(4) If red sites do percolate, the diffusion over the red 
cluster is possible and does not need activation. The 
temperature dependence saturates. 

(5) If red clusters do not percolate, the diffusion has 
to go via white sites, and therefore involves an activation 
process; its temperature dependence shows the Arrhenius 
behavior. 

These features, predicted by EMA, have also been 
found in simulations of a genuine ternary lattice in which 
red clusters run clung on the black chains by construc- 
tion. Figure [7| shows the behavior of Arrhenius plots for 
low polymer concentrations in the case of attractive in- 
teraction. It shows the logarithm of the normalized effec- 
tive diffusivity for different values of 03 . For 03 < 0.06, 
polymers remain sparse and isolated, their red perimeter 
sites don't percolate, no infinite red cluster exists and the 
system is in an Arrhenius regime (5). When the number 
of chains is increased, the transition from the Arrhenius 
to the saturation behavior (4) is observed at the critical 
value 03 = 0.06, revealing the emergence of an infinite 
red cluster. This critical value is far below the usual per- 
colation threshold of a cubic lattice due to the fact that 
red sites are arranged in connected groups on the perime- 
ters of black chains. This number can not be predicted 
by simple EMA and can be translated into an estimate 
of the percolation threshold of perimeter sites of chains. 



CONCLUSIONS 

We have considered diffusion of small molecules in a 
solid polymeric medium taking into account the inter- 
action between polymers and diffusing particles which 
can be both attractive or repulsive. The diffusivity has 
been analyzed from different perspectives both analyti- 
cally, using a modified effective medium approximation, 
and numerically by performing direct Monte Carlo simu- 
lations. While the diffusivity is only slightly affected by 
the chain's length, its temperature dependence crucially 




FIG. 7: (color online) Simulation graph: Arrhenius plots of 
the normalized diffusivity in the attractive case vs. e for dif- 
ferent segment concentrations. 



depends on the kind of interaction. This behavior de- 
pends on the sign of the interaction energy and is related 
to the existence of a percolating cluster of interaction 
sites surrounding polymer segments and/or a percolat- 
ing cluster of voids on which particles are free to travel 
without activation. 
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